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SUMMARY 

A small-scale laboratory magnetic suspension system, the Large Angle Magnetic 
Suspension Test Fixture (LAMSTF) has been constructed at NASA Langley Research 
Center. This paper first presents some recent developments in the mathematical 
modelling of the system, particularly in the area of eddy current effects. It is shown 
that these effects are significant, but may be amenable to modelling and measurement. 
Next, a theoretical framework is presented, together with a comparison of computed 
and experimental data. Finally, some control aspects are discussed, together with 
illustration that the major design objective of LAMSTF - a controlled 360° rotation 
about the vertical axis, has been accomplished. 


INTRODUCTION 

In order to explore and develop technology required for the magnetic suspension 
of objects over large ranges of orientation, a small-scale laboratory development system, 
the Large Angle Magnetic Suspension Test Fixture (LAMSTF) has been constructed at 
NASA Langley Research Center. Possible applications for magnetic suspension systems 
of this general class include space payload pointing and manipulation, microgravity 
vibration isolation and wind tunnel model suspension [1]. An important objective of 
this particular project is to investigate the dynamic modelling of large-gap magnetic 
suspension systems, so that future systems can be designed with higher confidence 
levels. 


HARDWARE DESCRIPTION 

The general configuration [2] is illustrated in Figures 1,2. An array of five, room 
temperature, copper electromagnets are equally spaced on a 13.77 cm radius. The coils 
are wound with 509 turns of AWG 10 enamelled copper wire on bakelite forms, with 
mild steel cores. The electromagnets are mounted on an aluminum plate 1.27 cm thick. 
Each electromagnet is driven by a transistor switching power amplifier, rated at ± 150 
V and ±30 A continuous, with four-quadrant operation. 

The suspended element consists of 16 wafers of Neodymium- Iron-Boron 
permanent magnet material, each 0.851 cm in diameter and 0.3175 cm thick, epoxied 
into an aluminum tube, 5.32 cm long and 0.9525 cm outside diameter. The total mass 
of the suspended element is 22.5 grams and the moment of inertia about tranverse axes 
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is 5.5 xlO' 6 kg.m". The direction of magnetization is along the axis of the cylinder, 
which is horizontal when suspended. The nominal magnetization is 954,930 A/m (1-2 
Tesla), although measurements have indicated a slightly lower working value The 
suspension height is 0.1m, measured from the axis of the suspended element to the top 
plane of the electromagnet conductor. 

The position sensing system consists of multiple light beams, arranged in the 
vertical and horizontal planes, partially interrupted by the suspended element The 
light sources are miniature infra-red light-emitting diodes, intended for use with fiber- 
optics, with collimating lenses added. The light receivers are matching infra-red 
phototransistors, with focussing lenses added. The complete sensor system is mounted 
on a framework which can be rotated by hand about a vertical axis. A schematic 
diagram of the sensor assembly is shown in Figure 3. 

Several different control systems have been developed and demonstrated, 
including a simple analog version with phase-advance ( lead ) compensation [2], a f i r s t - 
generation digital controller, generated using the bilinear (Tustin’s) transformation [3,4], 
a decoupled PD controller [5], and LQR and LQG designs [4]. 


EDDY CURRENT EFFECTS IN LAMSTF 
Introduction 


Whenever a time-varying magnetic flux penetrates a conducting medium, eddy- 
currents will be generated. In the case of LAMSTF, the principal time variation is due 
to the necessary control forces and torques being generated by fluctuating electromagnet 
currents, since the system is open-loop unstable. In the original design, eddy-current 
circuits were deliberately introduced in three main areas, as illustrated in Figure 4 . 

1) Position sensor structure, 2) Electromagnet cores, 3) Aluminum baseplate 

This was done so that it would be necessary to measure, analyze and model the 
eddy current effects, rather than attempting to avoid their influence, as is the usual 
practice. The fact that stable suspension was initially achieved rather easily [2J was 
taken to indicate that the eddy current effects were not very significant. However, a 
consistent discrepancy has been found in the dynamic behaviour in the pitch degree- 
of-freedom, illustrated in Figure 5 [3,4]. In consequence, an analysis and modelling 
effort has now been undertaken. 


Eddy Current Modelling 


A simplified analysis can be employed to assess the effects of eddy currents in 
LAMSTF. The two key assumptions are some a priori knowledge of the geometry of 
the eddy current circuit and that the circuit geometry be independent of frequency. 
The first assumption might require that the eddy currents be constrained to flow 
around well-defined paths, such as the position sensor structure here, rather than 
through large plates or shells of conducting material. Alternatively, the circuit 
geometry must be relatively simple and predictable. The second assumption requires 
that the “skin depth” be much greater than the local material thickness. The skin 
depth is given by the following formula [6] : 

d = y/ ( go Wv fna ) or ^ {ww*) ' (1) 
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- where d = Skin depth, p = Permeability, p = l/a — Resistivity, uo — angular 
frequency. In the case of LAMSTF, the natural frequencies of the suspended element 
are rather low, of the order of 10Hz or less. For an aluminum conductor, the value at 
10Hz would be around 28mm, much greater than the typical material thickness in 
LAMSTF. The only exceptions are the iron electromagnet cores, although it is found 
later that their influence is confined to higher frequencies. 

If both of the above assumptions are satisfied, the resulting model corresponds to 
that commonly described in literature as the Single Time Constant Model. The 
derivation resembles the analysis of a transformer with a shorted secondary, as 
illustrated in Figure 6 : 
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- where R e , L e are the resistance and inductance of the n’th eddy current circuit and 
L m is the 11 mutual inductance between the primary (the electromagnet coil) and the 
eddy current circuit. Note that mutual inductances between multiple eddy current 
circuits are neglected. The terminal characteristics of the primary (driven coil) can be 
found by combining these two equations : 
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One special case is of interest here. Suppose that : 

L = aL Cl ( 0<a<oo ) and L m =^(LL e ) ( 0 < /3 < 1) 
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This indicates that a non-dissipative (reactive) secondary effectively “turns off’ 
part of the primary inductor. Continuing, the field components generated (at the 
suspended object) can be expressed as : 


B i= K j I + V*. + ■ 
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where K-, K e are constants representing the field generated at the suspension location 
by the electromagnet and the n’th eddy current respectively. Now the factor K e will, in 
general, be different for each field component, that is each individual eddy current will 
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affect each field component by a different proportion. Therefore the eddy current 
effects in a system involving several electromagnets and eddy current circuits should be 
represented as follows : 


[Bj] = [K.) [I] + [KjJ [I e ] - (8) 

- where [B-] = (B x By B z B xx • •• )^, [I]=(Ij ^2 ' ) > [K-] is a rectangular matrix of 

J J 

field coefficients and [K ; ] is similar, though possibly of differing dimension. It is 

J 6 

presumed that [I e ] can be derived from [I], following equations 2,3. 

Alternatively, if the eddy current circuit has similar geometry to the primary 
(for example the induced current in electromagnet cores), it can be argued that the 
relative effect on all field and field gradient components at the suspended object will be 
similar. In this case, the representation can be considerably simplified by invoking a 
false current as shown : 



Kej L m ^ s 

Kj ( R e j + LejS ) 
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where Bj = Kjf 
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It should be noted that the change in electromagnet terminal characteristics and 
the change in field at the model are two separate effects and should be modelled as 
such. 


Determination of Parameters 


The question now is, can the parameters K e , L e , R^ and Ljp n be estimated 
and/or measured with sufficient accuracy? First the problem of estimation is addressed. 

Calculations have been carried out using the finite element computer code 
VF/GFUN, by Vector Fields Inc.. It should be noted that this code is magnetostatic 
and has no capability for direct eddy current calculations, although such codes are 
available (for instance ELEKTRA, by the same supplier). Instead, the code is used to 
calculate flux linkages, hence inductances, using : 


*n = | L ii 'j 


-( 10 ) 


VF/GFUN calculates the field on a grid representing the linkage plane of the 
eddy current circuit. The field normal to the plane is then numerically integrated (by 
the OPERA pre- and post-processor) to yield the flux linkage term. Figure 7 illustrates 
the general arrangement. The calculation of the K Cn terms is straightforward. 

By way of example, a series of calculations has been made for a single LAMSTF 
electromagnet with a representation of one part of the position sensor assembly 
mounted on the same axis, as illustrated in Figure 8. The required parameters were 
predicted (or previously measured) to be : 


L = 0.0275 H 
L e = 6.69 x 10' 7 H 
L m = 1.0998 xlO' 5 H 


R = 0.74 fi 

R e = 2.243 x 10' 4 n 
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K z = 3.495 x 10' 4 T K Ze = 4.369 x 10' 6 T 


Incorporating these values in equation 7, and examining the axial (z-axis) field 
component, gives : 


B z 



6.13 x 10~ 4 s \ 
1 + 2.983 x 10' 3 s ' 


-( 11 ) 


It is seen that the reasonant frequency of this eddy current circuit is around 
53Hz, significantly higher than LAMSTF open-loop natural frequencies, but still well 
within the range of interest. 


Experimental Verification 


Actual measurements of the current to field transfer function, corresponding to 
equation 7, were made with an experimental set-up as described above, and later with 
LAMSTF. Field components were measured with a F.W. Bell Model 9903 Hall-effect 
gaussmeter. Electromagnet currents were measured using a current shunt. The transfer 
function was measured directly with a Schlumberger Model SI 1250 analyzer, with sine- 
sweep excitation. The results for an air-cored electromagnet with no eddy current 
circuits are shown in Figure 9, and are taken to represent the probe -f- instrument + data 
acquisition system response. These results are subtracted from all subsequent 
measurements. Figure 10 shows the measured transfer function for B z , together with 
the predictions from equation 11. The agreement is thought to be satisfactory. The 
values of most parameters could be adjusted (refined) to give a better agreement, as 
shown in Figure 11. The only significant residual discrepancies are seen to occur at 
higher frequencies where the validity of the Single Time Constant Model is 
questionable. 

The adjusted form of equation 11 is : 


B z = 



7.591 x 10' 4 s \ 
1 + 2.934 x 10' 3 s ' 


- (lib) 


More Complex Cases 


If the electromagnet is mounted on the aluminum plate, a second eddy current 
circuit is added; when the iron core is inserted, a third is added. Figure 12 shows the 
comparison between experimental and computed responses. Again, the agreement is 
fair, although capable of improvement by refinement of parameter estimates, also shown 
in Figure 12. Note that, even if refinement of parameters is undertaken, the model does 
not correctly predict the high frequency behaviour, particularly where the iron core is 
present. This is due to the iron core becoming highly dissipative at these frequencies, 
due to its small skin depth. 

An additional series of calculations and measurements has been made for the 
vertical field component generated at the centroid of the suspended element due to a 
LAMSTF electromagnet at the design location, i.e. with the sensor ring off-axis relative 
to the electromagnet. The geometry is illustrated in Figure 13 and typical computed 
and measured responses are shown in Figure 14. It is clear that the model is 
progressively less suitable for progressively more complex geometries. 
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Measurements have also been made with the full LAMSTF sensor frame in place. 
Examples are shown in Figure 15. Due to the geometrical complexity, with multiple 
interlocking eddy current paths, it has not yet been possible to compute resonable 
estimates by the methods shown. It should be noted, however, that another assumption 
inherent in the form of the model previously chosen has been violated, that is the lack 
of interaction between separate eddy circuits. If these circuits are physically and 
electrically connected, this is clearly not reasonable. 


Terminal Characteristics 


It appears to be possible to experimentally estimate certain important 
parameters without direct measurement of magnetic fields. Figure 16 shows a 
comparison of measured and computed terminal characteristics for the single LAMSTF 
electromagnet mentioned above. The agreement is not perfect, but sufficient to 
validate the approach and can, of course, be improved by adjustment of parameters. 


PART 2 - IMPLEMENTATION of LARGE ANGULAR ROTATION 

Linear Modelling 


A portion of the development of the idealized governing equations for LAMSTF 
will be briefly reviewed here. Further details are given in References 7,8. The magnetic 
forces and moments acting on a magnetized core in a “large-gap” magnetic suspension 
system can be approximated as follows : 

F c ~V(M.VBo) T c ~V(MxB 0 ) (12) 

- where M represents the magnetization of the magnetic core in A/m, B the applied 
magnetic field in Tesla, V is the volume of the magnetic core in m , and the subscript o 
indicates that the field or field gradient is evaluated at the centroid of the magnetic 
core. Now, following the detailed development presented in reference 9, the effect of 
changes in relative orientation between the magnetic core and the electromagnet array 
can be incorporated as follows : 

F C = V [TnJMTnJ- 1 a (13) 

T C = V Sx ( [T m ] B) (14) 


Where a bar over a variable indicates magnetic core coordinates, [<9B] is a matrix of 
field gradients and [T m ] is the coordinate transformation matrix from electromagnet 
coordinates to suspended element (magnetic core) coordinates. The axis systems and 
some nomenclature are clarified in Figure 17. 


The fields and field gradients are created by an array of n electromagnets. Thus 
we can write : _ f 

B = ( K b ] r J — (15) 


max 


- where I = (L, I 2 , ••• , I n )^ and [ Kg ] represents a matrix of field coefficients. 
Similarly, each element of [<9B] can be written : 
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Iv 
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^®ij J Imax 


(16) 


At this stage, the equations are quite general, but will now be considerably simplified 
for the LAMSTF application. First, the magnetization of the suspended element is 
assumed to be along its principal (x) axis : 


M = ( Mjj ,0,0) 


(17) 


Continuing, the only large rotation of the suspended element is assumed to take place 
about the z axis, so [T m ] becomes : 


[T; 


mj 


Cos 9 Z Sin 6 Z 
— Sin 0 Z Cos 9 Z 
0 0 


Substituting, expanding and collecting terms : 


- where : 


T- 

± z 


Fz 


[KT] 

[KF] 


; kt 


0 

0 
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-LK Z J 

Sin 9 Z [ K x J + Cos 9 Z [ K y J 


(18) 


(19) 


(20) 


[ kf ] = 


Cos 2 9 Z [ K xx J + 2 Cos 9 Z Sin 9 Z [ K xy J + Sin 2 9 Z [ K yy J 
Cos 9 Z Sin 9 Z L K xx J + (Cos 2 9 Z - Sin 2 9 Z ) [ K xy J + Cos 9 Z Sin 9 Z [ K yy J 
Cos 9 Z [ K xz J -f Sin 9 Z [ K yz J 

- (21) 


Now there are two related problems to solve. Firstly, the equilibrium currents required 
to support the weight of the suspended element may be found by using : 


- where m c 


F Z = F z = m c g = V Mjj. 

- (12) 

T- = T- = F- - F- = 0 
x y L z r x r y u 

is the mass of the suspended element. Equation 19 can now be solved for 
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required currents as a function of yaw orientation 0 Z . As a point of interest : 


= 


m c g 


J xz - VMjj 


= 0.0962 Tesla/m in this case 


( 13 ) 


Secondly, a decoupling matrix is required so that the conroller can call for the 
appropriate combination of electromagnet currents to create the required forces and 
moments. This can be found by direct inversion of the [KT/KF] matrix from Equation 
19, giving : 


I 


demand 


[I<T] 

[KF] 


-1 


T- 

x z 


demand 

Note that the [KT/KF] matrix is a function of 0 Z . Initially, it has been found 
convenient to normalize the inverted [KT/KF] matrix column-by-column, with the 
normalization factor being incorporated into the loop gain. 

The variable decoupling matrix is incorporated into the controller in the 
following way. Matrices are calculated at frequent intervals (6° presently) and stored in 
memory. The controller interpolates in real-time between these matrices. The first 
three matrices, individually normalized, are shown here for reference : 


I 


demand 


0 = 0 ° 


I 


demand 


0 = 6 ° 



0.7713 

0 

1 

0 

-1 


1 

H 

1 


1 

1 

-0.8091 

0.6182 

-0.3092 


Tz 

= 

0.8584 

0.6178 

0.3091 

-1 

0.8092 


Fx 


0.8584 

-0.6178 

0.3091 

1 

0.8092 


F y 


1 

-1 

-0.8091 

-0.6182 

-0.3092 


Fz 


0.7637 

-0.1069 

1 

-0.2091 

-1 


T- 

y 


0.9695 

0.9340 

-0.6841 

0.7474 

-0.409 


Tz 

= 

0.8728 

0.6841 

0.1068 

-1 

0.7472 


Fx 


0.8234 

-0.5112 

0.5112 

0.8708 

0.8708 


F y 


1 

-1 

-0.9339 

-0.4091 

-0.209 


Fz 


demand 


demand 


374 


I 


demand 


0 = 12 ° 


0.7655 -0.2091 0.9339 -0.409 -1 


T- 

y 

0.9405 0.8708 -0.5111 0.8708 -0.5112 


Tz 

0.8916 0.7472 -0.1068 -1 0.6841 


Fx 

0.7958 -0.409 0.6840 0.7472 0.9339 


F y 

1 -1 -1 -0.209 -0.1069 


Fz 


demand 


For this interpolation, the controller must have information as to the actual yaw 
orientation of the suspended element. It was realised that under certain circumstances 
this information could be deduced in near real-time from the observed behaviour of the 
suspended element, rather than by direct global measurement of the suspended element 
orientation, in at least two ways. 


A. Yaw Error Tracking 

Suspension is established at a known reference orientation. If the sensor 
framework is now rotated (by hand) through some small angle, the suspended element 
will tend to become misaligned relative to the sensors. This yaw error signal is filtered 
to remove the effects of noise and the natural motion of the suspended element in 
response to small disturbances. Once the filtered error reaches a preset threshold, the 
controller steps through the array of decoupling matrices in an attempt to drive the yaw 
error back to zero. This method is rather crude, but has proven to be unexpectedly 
reliable. 


B. Current Distribution Analysis 

The predicted current distributions required to suspend at various yaw 
orientations, together with actual measurements, are shown in Figure 18. It should be 
noted that the distributions are almost perfectly sinusoidal in nature. By 
straightforward analysis of an observed current distribution, again with some filtering to 
remove the effects of noise and so forth, it is possible to deduce the orientation of the 
suspended element. 

In both cases, the algorithm cannot accommodate a steady-state applied yaw 
torque. Nevertheless, the fact that the position sensor framework does not require any 
form of angular orientation sensor is a great practical advantage. A sequence of 
photographs showing the system in operation is shown as Figure 19. 


DISCUSSION 

The simple eddy current model proposed appears to be satisfactory in the case of 
large eddy current circuits in conducting, non-magnetic material. Relatively simple 
computations are capable of providing reasonable estimates of important parameters, 
with the option of refinement based on measurements of magnetic field or 
electromagnet terminal characteristics. In the case of the iron electromagnet cores, or 
at higher frequencies, more elaborate models have been proposed (for instance [9]), but 
these have one potentially serious drawback. This is that the greater the elaboration in 
the eddy current model, the more complex the overall suspension system model 
becomes, and the greater the potential difficulty in manipulating this model in the 
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process of controller synthesis using modern model-based design methods. In some 
applications, the simple model proposed, with parameter adjustment, may adequately 
describe the eddy current influence on the dynamic behaviour, hence control 
performance, of the system. Discrepancies at frequencies well outside (above) the 
controller bandwidth would be of no consequence. 


CONCLUSIONS 

A simple model for the effect of eddy currents in the metallic structure of 
LAMSTF has been proposed and validated by experiment. However, it has not yet 
been possible to fully describe the eddy currents in the position sensor framework due to 
the geometric complexity involved. 

Eddy currents have been shown to seriously affect field and field gradient 
components in the frequency range of interest, such that they must be incorporated into 
a system dynamic model if modern control synthesis techniques are to be fully 
successful. 

The principal design objective of LAMSTF, the 360° rotation about the vertical 
axis, has been achieved. 
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Figure 7 - Illustration of OPERA Flux Linkage Computation (B z shown) 
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Figure 8 - Simplified Eddy Current Test Set-Up 
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Figure 14 - Typical Field at S uspended Element with Full (Aluminum) Sensor Frame 


Electromagnet Terminal Characteristics 



Electromagnet Terminal Char act eristics 



Figu re 15 - Electromagnet Terminal Characteristics 
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Practical Demonstration of Large Angular Rotation (0°) 
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Practical Demonstration of Large Angular Rotation (60°) 



Figure 18 (concluded) - Practical Demonstration of Lar ge Angular Rotation (12 0) 
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